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Abstract. The Landen transformation (a, b) ^ ((a + fe)/2, yab) preserves the 
value of an eUiptic integral and its iteration produces the classical arithmetic- 
geometric mean AGM(a, b). We present analogous transformations for rational 
functions integrated over the whole real line. 



1. Introduction 

The problem of indefinite integration of rational functions R{x) = B{x)/A{x) 
was finished by J. Bernoulli in the eighteenth century. He completed the original 
attempt by Leibniz of a general partial decomposition of R{x). The result is that a 
primitive of a rational function is always elementary: it consists of a new rational 
function (its rational part) and the logarithm of a second rational function (its 
transcendental part). 

In the middle of the nineteenth century Hermite [?] and Ostrogradsky [?] de- 
veloped algorithms to compute the rational part of the primitive of R{x) without 
factoring A{x). More recently Horowitz [?] rediscovered this method and discussed 
its complexity. The problem of computing the transcendental part of the primitive 
was finally solved by Lazard and Rioboo [?], Rothstein [?] and Trager [?]. For 
detailed descriptions and proofs of these algorithms the reader is referred to [?] and 
[?]■ 

This paper contains a method of computing definite rational integrals that, unlike 
the methods described above, does not involve the factorization of any polynomial. 
In this new method, the value of the integral is obtained as the limit of a sequence of 
transformations of the coefficients of the integrand. Thus, the algorithm presented 
here is in the spirit of the classical Landen transformation for elliptic integrals. 
These are integrals of the form 

dx 



(1.1) K{k) , , 

' ^(l-a;2)(l-fc2a;2) 

that have been studied since the eighteenth century. The reader will find in [?] 
more information about them. Its trigonometric version, 

G[a,b) = / , „ , 

JQ V 0^ cos2 9 + sin 9 

was considered by Gauss [?] in his work on the lemniscate. The special case k — i, 

dx 



(1.2) Lem = 



vT 
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appears as an expression for its arc length. He inferred from a numerical evaluation 
that the function G(a, b) is invariant under 

(1.3) €:(a,5) (i(a + 6), Vafe) . 

A transformation of the parameters of an integral is called a Landen transformation 
if it preserves the value of the integral. The example (jl.3p is the original one. 

It is a classical result that the iteration of £ produces two sequences, a„ and 
bn, that converge quadratically to a common limit: AGM(a, 6), the arithmetic- 
geometric mean of a and b. The invariance of the elliptic integral (jl.2p yields 

(1.4) G{a,b) 



2AGM(a,6)' 



Iteration of (|1.3p provides a method to evaluate the elliptic integral G{a,b). For 
instance, four steps starting at oq — 1, bo = \/2 yield 22 correct digits of the integral 
in p.2p . See [?] for details and its relation to modern evaluations of n. 
We consider here the space of rational functions 



:= i?(.)=^(") 



Aix) 



P p-2 

A{x) = J2 "■kx^"'' and B{x) = ^ ^fcX^^^-fe 

fe=0 fc=0 



We assume 

• The degree p is an even positive integer. 

• The coeffients and bk are real numbers. 

• The polynomial A(x) has no real zeros. 

Under these assumptions the integral 

/oo 
R{x) dx 
-oo 

is finite. 

We describe a transformation on the parameters 
(1.6) *Pp := {ao, ai,- ■ ■ ,ap; bo, bi,- ■ ■ ,bp_2} 

of i? G £Hp that preserves the integral /. In fact, we produce a family of maps, 
indexed by to e N, 

^m^p • ^ ^pi 

such that 

/oo />oo 
i?(a;)da; = / 2,n,p{R{x)) dx. 
-oo J ~OQ 

The maps £m,p induce a rational Landen transformation on the coefficients: 

(1.8) ^m,p ■■ ^ K^^- 

We provide numerical evidence that the iterates of this map converge to a limit, 
with convergence of order to. 

In the case to = p = 2, we will show that the integral 

f°° dx 

(1.9) 1(00,01,02) = / 5- ■ 

aox^ + aix + 02 
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is invariant under the transformation 

2aoa2 



(1.10) ao 



ai 



a2 



aa + a2 
Qi(Q2 - gp) 

ao + 02 
(ao + a2)^ - aj 



2(ao + a2) 

This example is discussed in detail in [?]. 

The theory of Landen transformations for rational integrands is divided into two 
cases, according to the domain of integration. 

Case 1: The interval of integration is not the whole real line. 

Integration over a finite interval [a, b] is transformed to the half-line [0, oo) by a 
bilinear transformation. In detail, 

f'' , . , . f°° /a + bt\ dt 

(1.11) / R{x)dx = (b-a) R' * 



1 + t ; (l + i)2- 

Similarly, integration over half-lines [a,oo) and {—oo,a\ can be reduced to [0,cxd) 
by translations and reflections. Thus, the interval [0,cxd) encompasses all integrals 
that fall in this case. 

Landen transformations for even rational functions on [0, oo) were established in 
[?]. For example, the integral 

n rrt h ^ ^ cx^ + d^ + e 

(1-12) U6[a,b;c,d,e) := / „r , , , -, dx 

is invariant under 

(1.13) a 



b 



afo -I- 5a -I- 56 -I- 9 

(a -f 6 -I- 2)4/3 

a + b + 6 
(a + & + 2)2/3' 



with similar rules for the coefficients c, d and e. 

The map (|1.13p can be iterated to produce a sequence (a„, 6„; c„, d„, e„) with 
the property 

(1.14) C/6(a„,6„;c„,d„,e„) = Ue{a,b; d, e). 

Its convergence was discussed in [?], assuming that the initial conditions ao, bo are 
nonnegative. The main result is the existence of a number L, depending on the 
initial data ag, • • • ,eo, such that a„ —^ 3, 6„ — > 3, and c„ — > i, d„ ^ 2L, and 
e„ — + L. The convergence is quadratic. 

The positivity condition on initial data was eliminated in [?] , where we reinter- 
pret the Landen transformation (|1.13p in geometric terms. The new integrand is 
the direct image of the original one under the map w = {z^ — \)/2z. In concrete 
terms, if R is the original integrand and 

(1.15) z±{w) = w± -J 'Up- 4- 1 
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are the two branches of the inverse of w, then the new integrand is given by 

(1.16) ij(^^(^))^+i?(^_(^))^. 

This geometric interpretation extends to the algorithm presented in [?], where an 
analogue of (|1.13[) is given for an arbitrary even function. These transformations 
on the coefficients define a map, 

(1.17) $2„ ^]R2n-l^ 

which is the rational analogue of 2; in (|1.3p . These are named even rational Landen. 
Using this approach, we have established a proof that the iterations of $2n converge 
precisely when the initial integral is finite. 

A purely dynamical proof of convergence of the iterations of $2n is presented 
in [?], but only for the case of degree 6. The relation between (|1.13p and the 
invariance of the rational integrals is still part of the argument. It is established 
that the iterations are eventually mapped to the first quadrant, and then the results 
of [?] are applied. It would be desirable to obtain a proof of convergence completely 
independent of the integrals that gave origin to these maps. 

The existence of this type of transformation for an odd rational integrand is an 
open question. 

Case 2: the domain of integration is the real line. 

This is the case we present here. We give a Landen transformation for integrals 
over M. The convergence of the iterations of these maps can be established along 
the lines of [?], but a more direct analysis is still an open question. The issue of 
convergence is not discussed here, except for the numerical examples in Section [7l 

The new integrands, 2,m.p{R{x)), depend on the parameter p, the degree of the 
denominator of the original integrand, and the parameter m, the order of conver- 
gence. Both parameters are arbitrary. 

Section 2 presents a preliminary example that illustrates the methods developed 
in the rest of the paper. Section 3 introduces two families of polynomials that 
are the basis of the rational Landen transformations. Section 4 consists of some 
simple trigonometrical identities. The integrand is scaled in Section 5, using the 
polynomials studied in Section 3. The algorithm leading to the rational Landen 
transformation is a consequence of the vanishing of a class of integrals. This is 
presented in Section 6. Examples are given in the last section. 



2. An example 

We begin with an example of a Landen transformation that introduces the meth- 
ods described in later sections. 

The integral of the rational function 



(2.1) R{x) 



x^ + 6x3 _,_ 29a;2 -|- 60x + 100 
is evaluated as 

38 TT 



/OO 
R{x) dx = 
-OO 



31a/3T' 
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using the factorization 

(2.3) + 6a;^ + 29^2 + 60a: + 100 = {x^ + 3a; + 10)^. 

We will produce a new rational function, 

, . , , ^ 202x^ + 45a: + 97 

( . ) L2Amx)) - 400a:4 + 1080a;3 + 2969a;2 + 3024x + 3136' 

and show that it satisfies 

/OC f-OO 
£2A{R{x))dx = / R{x)dx. 
-OO J —OQ 

(The notation £2,4 indicates the degrees of the transformation used to produce this 
new function. Details are given in Section [6]). 

The first step is to multiply the denominator, 

(2.6) A{x) = a:-* + 6a;3 + 29a;2 + 60a: + 100, 

by 

(2.7) Z{x) = 1600a;^ - 960a;^ + 464a:2 - 96a: + 16, 

so that E{x) ~ A{x)Z{x) can be written as a homogeneous polynomial in the 
variables 

(2.8) P2{x)=x^-l and Q2(a:) = 2a;. 

(These polynomials will be described in Section O) In detail, 

4 

(2.9) E{x) = ^e,P2'"'(x)Q^(a:), 

with eo = 1600, ei = 4320, 62 = 11876, eg = 12096, and 64 = 12544. Then, with 
C{x) — B{x)Z{x), we obtain 

(2.10) / = r^dx=r^dx. 



A{x) E{x) 

Now write 

(2.11) E{x) = Qt{x)[Y,eiR2^ 



vi=0 



where 

^o^o\ ( \ -P'2(x) x^ - 1 

(2-12) -^2(2;) = 77-^ = . 

Q2[x) 2x 

We would like to make the change of variables y = R2 (x) in (|2.10p . The function 
R2{x) has a multivalued inverse, with its two branches given by 

(2.13) X = y± Vy' + l. 

Therefore, we must split the evaluation of the original integral at the singularity 
a; = of R2{x). The identity (|2.10p is written as 

r C{x) ^ p C{x) ^ 



Ei{y) J -00 Ei{y) 
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where 

4 

(2.14) E,{y) = Y.^iy^~' 

1=0 

and 



(2.15) N,iy) = ^I't'f^. ^iy^^^)- 

The new integrand, {N+{y)+N_{y))/ Ei{y), corresponds to the expression in (|1.16p . 
A direct calculation shows that 

(2.16) N_{y) + N+{y) = 4(202y2 + 45?/ + 97), 
so that 

.2l7^ J ^ r 202^2 + A5y + 97 

^ ■ ^ 400y4 + 1080?/3 + 2969?/2 + 3024y + 3136 ^ 

202y2 + 45?/ + 97 
_^ (20y2 + 27y + 56)2 ^' 

as claimed. 

A proof of a transformation of this type for a general rational integrand is pro- 
vided in the next four sections. 

3. A FAMILY OF POLYNOMIALS 

For m e N, we introduce the polynomials 

Lm/2J 

(3.1) PM = E (-i)'(^j^""'' 



J=0 
L(m-1)/2J 

a; 



2j + l 



jl ^ \ m-(2i+l) 



(3.2) Q,„(a:) = E (-!)■'" 

which play a fundamental role in the algorithm discussed here. They will comprise 
the numerators and denominators of a natural change of variables discussed in the 
last two sections. 

The degrees of Pm and (5,„ are m and to — 1, respectively. Observe that 

(3.3) P2{x)=x^-l and Q2{x) = 2x 
have appeared in Section [51 

Proposition 3.1. Let M{x) = and fm(x) — x"\ Then the rational function 
Rm = Pm/Qm Satisfies 

(3.4) Rra - M-'of^oM. 
Proof. The identity follows from 

{x + i)"" + {x ^ i)'"' ^ 2Pm{x) and {x + i)"" - {x - i)"' = 2iQ^{x). 

□ 
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Corollary 3.2. The function Rm satisfies 

(3.5) i?™ (cot 61) = cot(m6i). 

Proof. Use M(cot6') = e^''' in □ 

Note 3.3. The multiplicative property i?„ o = Rnm shows that the functions 
Rjn form a family of commuting rational functions. The cotangent function in 
appears as the limiting case of the Weierstrass elliptic p-function, 

1 1 



(3.6) pix) - ^+ J2 



{x - riiLUi - n2LU2)^ {niUJi + n2LU2) 



where the term ni = n2 = is excluded from the sum. In the case uji — 1 and 
UJ2 — > oo, we get 

d TT^ 

(3.7) p(x) ^ -vr— cot(^x) - —. 

The function p(nx) is even and elliptic, therefore it is a rational function gn of p. 
In view of g^ o g^ = gnrm these functions commute. An extraordinary fact, due to 
Ritt [?], is that these are all such commuting rational maps. The functions _R„ are 
a special class of the gn- See [?], section 2.13, for details. 

The identity (|3.4p permits the explicit evaluation of the zeros of Pm and Qm- 

Proposition 3.4. The polynomials Pm and Qm have simple real zeros. Those of 
Pm are given by 

/(2fc+l)7r\ 

Pk = cot -^^ — for < fc < m - 1, 

\ 2m / 

and those of Qm are 

f k'n\ 

qk — cot ( — ) for 1 < /c < TO — 1 . 
\m J 

Proof. The identity Rm — o fm o AI yields 

Rm{Qk) = M-V„(M(cot(fc7r/TO))) 

= M-\fmie"''''/n) 
= M-\l) = oo, 

so that Qm{gk) = 0. The degree of Qm is to — 1 and the qk are all distinct, hence 
these are all the zeros. The argument for pk is similar. □ 

The polynomials 

Lm/2J 

(3.8) P;.(a) = J2 (-1)1 J ^^''^'^'^ 



i=0 
L(m-1)/2J 



(3.9) QM = J2 (-l)H2z + l 

have appeared in our development of definite integrals related to the Hurwitz zeta 
function. See [?] for details. They are connected to Pm and Qm via 

(3.10) P;,(a;) =x"P™(a;-i) and Q*m{x) ^ x"'-^Qm{x-'). 
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The role of these polynomials in the development of the Landen transformation 
comes from their trigonometric properties. 

Proposition 3.5. The polynomials P,* and Q*„ satisfy 

/ ^^ cosmd , ^» / „^ smm6 

3.11 P* tan0 =^— and Q;^{tan9) ^ 

cos'" cos™ 

Proof. We give the details for Q^. The series expansion 



, . smjx tan ^ t) _ ^ (~1)^ (x)2fc+i 2fc+i 

^ ^ ^ (l+<2)x/2 (2fc+l)! 



fc=0 

where (a;)^ = a;(a;+ l)(a; + 2) • • • (x + k—l) is the Pochhammer symbol, is established 
by checking that both sides satisfy the equation 

(3.13) (i + f2)A^2i(x + l)^+x(a; + l)5 = 0, 
with the initial conditions .g(0) = 0, .g'(0) = x. Then (a;)*; reduces to 

(3.14) (-m)„ - 
for n < m, and vanishes for n > m, since m is an integer. Thus p.l2p reduces to 

L(m-1)/2J 

,2fc + 1, 



sin(mtan-^t) = (l + ^-Z^ ^ ("1)' JIJ^ 



fe=0 



2k+l 



for X — —m. This is equivalent to the second formula in (|3.1ip . A similar argument 
establishes the expression for P^. □ 

In terms of the original polynomials, (|3.1ip becomes 

„ , cos m9 , „ , sinmO 

3.15 ^'m cot^ and Q„ cot0 = ^— . 

sm sm 



4. A TRIGONOMETRIC REDUCTION 

The example described in Section [2] can be extended by using the transformation 
y = Rm{x) with higher values of m. The explicit evaluation of the new integrals 
requires knowledge of the branches of the inverse map x = R^{y). This is imprac- 
tical for m > 3. An alternative method is described in the next section. 

The explicit formula for the Landen transformation uses an expression of sin" cos^ 0, 
for a, 6 £ N, as a linear combination of trigonometric functions of multiple angles. 

We introduce the notation 

(4.1) c= l^'] andd= LfJ- 

The reduction formulas given below are expressed in terms of the function 

Some of the identities presented here can be found in the table appearing in [?], 
page 30. 
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Proposition 4.1. Let a, b G N and u e M. Then sin" u cos'' u is given by 
(-1)'^ 



2a+b 

HI! 

2a+b 

2a+h 

(-1)" 
2a+b 

Proof. Start with 
(e^"-e-™)''(e' 



T,(a, + ^ {Tc+j{a, b) + T,_j(a, 6)) cos(2ju) 



i=i 



^ (re_i+,(a, 6) + r,_,(a, 6)) cos((2j - 1) 



^ (re_i+,(a, 6) - r,„,(a, 6)) sin((2j - l)^.) 



X! (^c+j(a, fe) - Tc-j{a, b)) sin(2ju) 



for a even and 6 even, 



for a even and 6 odd. 



for a odd and b even. 



for a odd and b odd. 



a b 



iu{2k—a) 



E 



Ju{2j-b) 



Therefore 



■ — n a-\-b a+6 



(4.3) sin'^.cos^u = E E(-ir ''^M " 



fe=0 1^=0 



V — k j \ k 



^iu{2u—a—b) 



The result follows now by eliminating the imaginary terms on the right hand side 
of (ESI). □ 



5. The scaling of the integrand 

In this section we describe a construction of the polynomials Z{x) and E{x), 
introduced in Section O These are used to produce an appropriate scaling of the 
integrand in 

B{x) 



(5.1) 



/ = 



A{x) 



dx, 



so that the new denominator is E(x). Recall that E is the homogeneous polynomial 
in the variables {Pm{x),Qm{x)) introduced in Section [H 

We now express the coefficients of E and Z in terms of those of A. This requires 
the elementary symmetric functions 

of the p variables yi, ■ ■ ■ ,yp. These are defined by the identity 



(5.2) 



\{{y-yi) = Y.i-iyo\^\y,,--- ,y,)y 



p-i 



1=1 



1=0 
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(5.3) A{x) = ^flfe^C- 



Theorem 5.1. Let p, m G N, and 

p 

k=Q 

be a polynomial with real coefficients. Then there exist mp + 1 coefficients, 
(5.4) 

with r = p{m — 1), such that 

(5.5) A{x)Z{x) = E{x), 
where 

r p 

(5.6) Z{x) ■.= Y,^kX^''' and ^(x) := [P„,(x)]P-' [Q,„(x)]' . 

k=0 1=0 

The coefficients ei are polynomials in the coefficients — , • • • , — . 

Note 5.2. The effect of the theorem is to scale the integrand B{x)/A{x) to C{x) /E{x), 
where C{x) = B{x)Z{x) and E(x) = A{x)Z{x). The degrees are recorded here: 

(5.7) deg(A) = p, deg(i3) — p~2, deg(Z) = r = pm — p, 

deg(C) = s = pm — 2, and deg(i?) = pm. 

Proof. Let {xi, X2, . . . , Xp} be the roots of A, each written according to its multi- 
plicity, so that 

p 

(5.8) A{x) = aoY]_ix ~ Xj). 

j=i 

The rational function 

(5.9) = ^ 

introduced in Proposition [331 is well-defined at all the roots Xj. This follows from 
the fact that the roots of Qm are real and our assumption that the roots xi, - ■ ■ ,Xp 
of A{x) — are not. For < I < p, define 

p 

(5.10) ei :=aS'(-l)'n^™(^j) ^ '^l^HRmixi), R^ix2) , . . ■ , Rmixp)) 
and the polynomial 

p 

(5.11) i?(a;) =^e,a;f-'. 

1=0 

We now consider the identity, 

(5.12) f[{y - R^{X,)) - f](-l)V(P) {Rra{xi), Rra{Xp)) v"'', 
3=1 1=0 

that comes from (j5.2l) . Clearing denominators, we obtain 

(5.13) n (Q„(a;,)2; - P„,(x,)) - '» e^yf-' = ac7™i?(y). 
i=i i=o 
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In particular, 

p p 

(5.14) H{Rm{x)) = a™ J] Q„,(x,) x J] - Rra{xj)) . 

Finally, define the polynomial 

E{x) = j2eiPP^-\x)Ql{x) 

1=0 

= H{Rr^(x))Ql,{x). 

Identity (|5.14[) shows that the zeros of E are precisely the values Rm{xj), I < j < P- 
The coefficients of E, given in (|5.13p . are symmetric polynomials of the roots Xj 
of A. The fundamental theorem of symmetric polynomials [?] states that e/ is a 
polynomial in ' ' i This, in turn, proves that e; e M and thus E G M.[x]. 

Now observe that (|5.14p yields E{xj) — and the corresponding factor Rm{x) — 
Rm(xj) appears with the same multiplicity as Xj. We conclude that A divides E 
and define Z to be the quotient. 

□ 



6. The reduction of the integrand 

In this section, we produce explicit formulas for rational Landen transformations 
of the integral 

, ^ r B(X) 

(6.1) / = / -TTT^^;. 



A{x) 

As before, we assume that A, B G M.[x], relatively prime, and that I < oo. 
In Section [5] we have scaled the integrand in (16. 1|) to the form 

(6.2) 7 = r^.dx, 
where the denominator is written as 

(6.3) E{x) = j2eiPl-Hx)QUx). 

1=0 

Here, P™ and Qm are the polynomials discussed in Section [S] The scaling of the 
denominator is achieved through multiplication by Z(x), as given in Theorem 15. II 
The numerator becomes 

s 

(6.4) C{x) = B{x)Z{x) ="^Ckx''-''. 

k=0 

The coefficients Ck are given by 

3 

(6.5) Cj = ^ Zkbj^k for < j < s, 

k=Q 

where 6i = 0ifi>p— 2 and Zi — if i > r — pm — p. 
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The change of variables x = cot 6 and the relations (|3.15|) yield 
where 

s 

(6.7) CT^^p{0) = ^ cfc cos^"'' 6* sin'^' 6* 
and 

p 

(6.8) ET„,p(0) = ^ e; cosP^' (to6I) sin' (m6l). 



1=0 



The discussion of this integral is divided according to the parity of m. Recall 
that p is assumed to be even. The details are presented in the case m odd. 
The parameter s — mp — 2 is even and we write s = 2X. Split (|6.6p as 



_ ^ r sin^J cos'-^J dd r sin^J+i 9 cos'^-^J-i Q dO 

and consider the evaluation of each of these integrals. 



The evaluation of Ii. The identity in Proposition 14. II yields 

sin2^^ cos^-2j" = I_il!r^(2j,s-2j) + 
2'' 



2^ 

k=l 

and replacing this in the definition of Ii yields 
(6.9) 

^ A. ^ ... r d9 

h — 



^ J2 iTx+k{2j, s - 2j) + Tx-ki2j, s - 2j)] cos(2fc0), 



lg(-l)^c,T.(2,,.-2,)^\^^^^^^^^^ , 

+ ^ E(-l)^-. E [^-(2.- ^ - 2.) + T._.(2,, . - 2,)] £ 

The periodicity of the integrand, and the fact that mp is even, show that the 
integrals appearing in (|6.9p are half of the corresponding ones over the whole period 
[0, 27r]. Thus, 



(6.10) 



1 , /-St 

A A „2ir 



1 V^. ^, V^r / ^ . M f '^'^ C0s(2ke) d9 

— Y.i^iyc2j J2 [T^A+fe(2j, s - 2j) + n-k{2j, s - 2j)] / ^ . 
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We now show that most of the integrals in ()6.10p vanish. This reduction is re- 
sponsible for the existence of a rational Landen transformation. 



Introduce the notation 

sin(fc6l) d0 



(6.11) S^,p{k) := 
and 

(6.12) CmAk) ■■= 



^Tm,p{0) 



cos{k9) de 
ET„j_p(6') 



Lemma 6.1. Let k, m, p G N be arbitrary. Then Sm.p{k) and Cm_p{k) vanish 
unless fc is a multiple of m. 

Proof. In the definition of Sm^p{k) let 9 i-^ 9 + 27rj/m for j = 0, 1, • • • , m — 1. The 
average of these m integrals is 

^ l'2iT ^0 I m-l m-l 

Sm,p{k) = — —= — sin(fc6') cos(27rfcj7m) + 008(^6*) sin(27rfcj7m) 

mJo ET„^p(6') \^ ^ ^ 

If k is not a multiple of m the integrand vanishes because the sums in it are the 
real and imaginary parts of 

™-l -I _ 27ri/c 

E^2Trijk/m _ ^ ^ _ q 

0=0 

A similar proof follows for Cm,p{k)- □ 



In view of Lemma |6. 11 we replace k by am, where 1 < a < — 1, with v — p/2. 
Then (pTU]) becomes 



(6.13) 



+ :^E(-1)^2,E[Ta+„™(2j,s-2j) + Ta_.™(2j,5-2j)] / ' ' 
The change of variables Lp = mO produces 

l'2-K 1 f-2'jTin f'2-K 

(6.14) / ■■■ de ^ — ■■■ d(p= ■■■dip, 

Jo rn Jo Jo 

using the periodicity of the integrand. We conclude that 
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(6.15) 

1 v^. N , . X . M f^"" cos(2a9) d0 

+ ^E(-l)^2,E[TA+»„(2j,5-2j) + rA-„™(2j,s-2j)]y^ ET,p{0) ' 



j=l a=l 

where the denominator is 



(6.16) ETi,p(6l) = ^e;cosP-'6'sin'( 



The next step is to bring back the domain of integration to [0, tt] . The symmetry 
of the integrand shows that the integral over [tt, 2tt] is the same as that over [0, tt]. 
We conclude that 



(6.17) 

1 x-^/ ^ ., ^ r„ , s „ / M co8(2a6) d6 

+ :^E(-1)^2,E[Ta+„™(2j,s-2j) + Ta-o™(2j,5-2j)] / . 

j=l a=l i^^l,p[0) 

The change of variables y = cot 6 gives, recalling that — p/2, 

Jo ETi,p(6l) H{y) 
where the polynomial 

(6.19) H{y) = Y.^iyP'' 



1=0 



was introduced in (|5.1ip . The identity p.lSp is now used to change variables in the 
second integral to obtain 



The next step is to write P2a{y) in terms of 1 + j/^. 
Lemma 6.2. The polynomial P2a{y) can be written as 

(6.21) P.„(y) = f](-l)"-^_^2M(^2«-A(i + ^2), 
[3=0 a p \ P J 
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Proof. Start with 

« /o„A 



3=0 ^ ■' 



X 



= E(-i)^(2-)[(i+-^)-i: 



a-j 



3=0 

I a~l3 



13=0 \ j=0 

and the result follows from 

Q-/3 



This sura arises as a corollary of Gauss's hypergeometric evaluation [?], 

V{c)T{c - a~h) 



(6.23) 2Fi[aMc\l] 



T{c-a)T{c-b) ' 



valid for Re (c — a — 6) > 0. In our case, a=-^ — a, 6 = /3 — a and c = ^, so that 
c~a — b — 2a~(i>Q. See [?], page 66 for a proof of (|6.23p . 

□ 

We return to the evaluation of Ji. The expression in ()6.17|) becomes 

h - ^|:(-i)'=«r.(2i,,-2i)£(i + !,^r-^ + 

A u-1 

+ ^ E(-l)'^2, E (^A+am(2j " S - 2j) + Tx-o^W^j " S - 2j)) X 



j=0 a=l 

,^0 2a-H P Jj-J'^'^ H{yy 



Expanding the powers of 1 + y^, and reversing the order of summation, leads to 

7=0 ^ ' ^ 3=0 •' ^' 

7=0 \j=0 a=l 13=0 J J-c 



y 



H{y) 

dy 
H{y) 



^E E E E Ml (j, p) / ; 

7=1 yi=0 Q=i'-7 /3=Q-i/+7+l / °° 



27 dy 

y ^ 



H{yy 
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where 

X [rA+am(2j, S - 2j) + rA_a,n(2j, S - 2j)] . 

The evaluation of l2- A similar calculation leads to 



^ /a-1 1^-1-7 a-l \ „c 

^2 = E E E^^2(j'"'/5'-'^') / 

+ ^E E E T.M2{j,a,(3;rn,p) \ / 

7=1 \j=0 a=iy—y fj=0 / 



" ^27+1 

oo^ H{y) 

" „27 + l^ 



where 

M2{j,a,f3;m,p) = (-l)^+'^C2,+i2 (2/3+1/ V 7 

X [Tx+c.m{2j + 1, S - 2j - 1) - Tx-cm{2j + 1, s - 2j - 1)] . 

For the convenience of the reader we summarize the information as a theorem. 



Theorem 6.3. Let p, to g N and assume p is even and to is odd. Define 

s p 

(6.24) s — mp — 2, r — p{m — 1), A = -, and ^ — 2' 

and consider the polynomials 

p p-2 



(6.25) A{x) = J2 akx^-'' and B{x) = ^ bkxP- 



-2-k 

'kx- 

fe=0 fc=0 

Then 



y-oo^(a^) J-ooH{x) 
The new denominator iJ is given by 

p 

(6.27) il(a;) = JI^'^""'' 

1=0 

where, for < j < p, the coefficients ej are solutions to the system (|5.5p . An 
expression for ej in terms of the coefficients ak is given in ()5.10p . 
The new numerator J employs the function 

X 

(6.28) T,ia,b) = ^(^1)--+^- 
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and is given by 



1/-1 / / i\ A 



(6.29) J(x) ^ i ME(-1)'^2,Ta(2j- .-2j) 

I/-2 / A 1^-1-7 a \ 

+ E E E^'^iO."'/5;m,p) x^-^ 

7=0 \j=0 c(=l /3=0 y 

+ ^EIE E E Mi(j,a,/3;m,p))x^'^ 

7 = 1 \ J=0 0=1/— 7 /3=a — 1/+7+1 / 
iy-2 /a-1 1^-1-7 Q-1 \ 

+ ^E E E E^^2(j'"'/5;-'^)h'''^^ 

7=0 yj=o Q=i /3=o y 

. iy-2 /a-1 u-1 q-1 N 

+ ^E E E E^^20'«'/?;™'P) 
7=1 yj=o Q=iy-7 0=0 y 



The coefficients Cj are given by 

3 



(6.30) Cj — E] ^kbj^k for < J < s, 



fe=0 



with bi=Oifi>p— 2 and = if i > r = — p. The values of Zj are obtained 
as solutions of the system (|5.5p . 



Finally, 

MiO,a,/3;™,p) = (-1)^+ ^ j ^ j 

X [rA+am(2j, S - 2j) + Tx-c.m{2j, S - 2j)] , 

and 



M2{j,a,P;m,p) = (-1)^+Pc2,+i2 



2/3+1, 

X [TA+a„(2j + 1, s - 2j - 1) - Tx-a.m{2i + 1, s - 2j - 1) 



Note 6.4. Surprisingly, the expressions for H and J given in (|6.27p and (|6.29p 
remain valid when m is even. This results from a similar calculation whose details 
are omitted here. 



7. Examples of rational Landen transformations 

This section contains some examples that illustrate the rational Landen trans- 
formations. 
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Example 7.1. We calculate the transformation for the case p = m = 2. The 
integrand in 

(7.1) / - / 2T ^i*^ ao^O, 

is quadratic, thus p = 2. We construct the Landen transformation with convergence 
order m = 2. Therefore, s = mp —2 = 2 and r = p{m — 1) = 2 in the notation 
defined in Theorem 15.11 The scaling of Section [5] amounts to finding parameters 
zq, zi, Z2, and eo, ei, 62 such that 

{aox'^ + aix + a2){zox'^ + zix + Z2) = eoPi{x) + eiP2{x)Q2{x) + e2Ql{x), 

with P2{x) = x^ — 1 and Q2{x) = 2x. The linear system (15. 5p is of order mp+ 1 = 5 
and we choose the free parameter 

4a2 

(7.2) eo = Q2{xi)Q2{x2) = AxiX2 = , 

ao 

according to the convention in (|5.10p . The solution of (|5.5p yields 

4a2 4ai 4 

(7.3) Zq = — , zi = ^2 = — , 

Oq Oq ao 

and the formulas in (|5.10p produce 

4a2 2ai(a2-ao) (00 + 02)^-0? 

(7.4) eo = , ei = ^ ' ^2 = 5 ■ 

ao Oq aQ 

Therefore, the denominator of the new integrand is 

(7.5) H{x) = 4a2^2^2ai(a2-ao)^^(ao + a2j^-af_ 

ao On Or, 



The new numerator is obtained from the formulas given in Theorem 16.31 In this 
case A = 1 and v = \, thus only one sum contributes to its value: 

(-7 R\ T/ \ 2(ao + a2) 

(7.6) J{x) = 2 • 

"0 

We conclude that 

dx f°° 2(ao + 02) dx 



_oo aox'^ + aix + a2 J ^aQa2x'^ + 2ai(a2 - ao)a; + [(ao + ^2)^ - a\] ' 

and (|7.ip is invariant under the transformation 

2aoa2 



(7.7) ao 

ai 
a2 

This was announced in (|1.10p . 



ao + a2 
ai(a2 - ao) 

ao + a2 
(ao + a2)^ - aj 

2(ao + a2) 
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Example 7.2. We present the Landen transformation of order 3 for the rational 

function 

(7.8) R{x) - ^,6 _^ iQ^5 + ii4a;4 _^ 4^^2x^ + 1041x2 + 1300x + 676 ' 

This is an example that violates the main assumption on the nature of the roots of 
A. Indeed, 

(7.9) A{x) = {x + 2f{x'^ + 6x + 13f and B{x) ^ {x + 2f, 

so that A has real roots. Although R is not reduced, it is integrable over R : 

y„oo (a;2 + 6a; + 13)2 ~ 16' 

This example shows that the rational Landen transformations preserve the existence 
of real poles of the integrand. Moreover, the real zeros that cancel these singularities 
are transformed accordingly to preserve convergence. 
The roots of A{x) — Q are 

(7.10) a;i = a;2 = —3 — 2i, X'i — x^^ —3 + 2i, x^ = xq = —2. 
The value of eo given in (|5.10p yields 

6 

(7.11) eo = Q3(^fe) = 269353744. 

k=l 

We have that p = 6 and m = 3, and so r = 12 and s = 16. Solving a system of 
order 19 yields 

(7.12) Ai{x) = (11a; + 2)2(373a;2 + 594a; + 481)2 
as the new denominator, and the new numerator is 

(7.13) Bi{x) = (11a; + 2)2(854x2 + 3240a; + 10709). 

Observe that the algorithm preserves the existence of a real root, but the root at 
X = —2/11 is cancelled. The reader will check the invariance: 



(7.14) / 

J — t 



854x2 + 3240x + 10709 , tt 

fix, 

(373x2 _^ 5942; _^ 431)2 iq 



Example 7.3. Finally we present a numerical example that illustrates the con- 
vergence of the iterative transformations constructed in this paper. The original 
integral is written in the form 

(7.15) / . ^ r + + + ■ ■ ■ + ^o'>-2 dx. 

flo j-co + flixP 1 + ao a2xP 2 + . . . + Qq ap 

The Landen transformation generates a sequence of coefficients, 

(7.16) {a("\a("',...,a(");6[,"\6i"\...,5("J2}, 
with *Po = as in (|1.6p . We wish to show that, as n ^ oo, 



(7.17) 



a2 "l "2 "p-2 

.ao a^ ' 6^ ' 6^ ' ' ^ 
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converges to 
(7.18) ( 0, 



.0, 



.0, 



9-1 



where q = p/2. The invariance of the integral then shows that 



9-1 
q-1 



(7.19) 



A) 

An) 



1 



The convergence of v := Un — Uoo to is measured in the L2— norm, 



(7.20) 



M2 



1 



/2p-2 

El 



1/2 



. fe=i 



V2p^ 
and also the Loo-norm, 

(7.21) = Max{|jt;fc|| :l<fc<2p-2}. 

The rational functions appearing as integrands have rational coefficients, so, as a 
measure of their complexity, we take the largest number of digits of these coeffi- 
cients. This appears in the column marked size. 

The following tables illustrate the iterates of rational Landen transformations of 
order 2, 3 and 4, applied to the example 

3a: + 5 



(7.22) 



/ = 



14x3 + 743,2 _^ 1842, _,_ 208 12 ' 



The first column gives the _L2-norm of u„ — Uoo, the second its Loo-norm, the third 
presents the relative error in (|7.19p . and in the last column we give the size of the 
rational integrand. At each step, we verify that the new rational function integrates 
to -77r/12. 

Method of order 2 



n 


L2-norm 


Loo-norm 


Error 


Size 


1 


58.7171 


69.1000 


1.02060 


5 


2 


7.444927 


9.64324 


1.04473 


10 


3 


4.04691 


5.36256 


0.945481 


18 


4 


1.81592 


2.41858 


1.15092 


41 


5 


0.360422 


0.411437 


0.262511 


82 


6 


0.0298892 


0.0249128 


0.0189903 


164 


7 


0.000256824 


0.000299728 


0.0000362352 


327 


8 


1.92454 X 10-* 


2.24568 X 10"* 


1.47053 X 10-* 


659 


9 


1.0823 X 10-16 


1.2609 X 10-16 


8.2207 X 10-" 


1318 



As expected we observe quadratic convergence in the L2— norm and also in the 
Loo— norm. The size of the integrand is doubled at each iteration. 
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Method of order 3 



n 


L2-iiorm 


Loo-norm 


Error 


Size 


1 


15.2207 




20.2945 


1.03511 




8 


2 


1.97988 




1.83067 


0.859941 




23 


3 


0.41100 




0.338358 


0.197044 




69 


4 


0.00842346 




0.00815475 


0.00597363 




208 


5 


5.05016 X 10- 


-8 


5.75969 X 10-* 


1.64059 X 10" 


-9 


626 


6 


1.09651 X 10" 


23 


1.02510 X 10-23 


3.86286 X 10- 


24 


1878 


7 


1.12238 X 10" 


70 


1.22843 X 10"™ 


8.59237 X 10" 


71 


5634 


Method of order 4 


n 


i2-iiorni 


ioo-norm 


Error 


Size 


1 


7.44927 




9.64324 


1.04473 




10 


2 


1.81592 




2.41858 


1.15092 




41 


3 


0.0298892 




0.0249128 


0.0189903 




164 


4 


1.92454 X 10" 


-8 


2.249128 X 10-** 


1.47053 X 10- 


-8 


659 


5 


3.40769 X 10- 


33 


3.96407 X 10-33 


2.56817 X 10" 


33 


2637 



8. Conclusions 



We have presented an algorithm for the evaluation of a rational integral over R. 
Numerical evidence of its convergence is presented. 

Acknowledgments. The work of the second author was partially funded by 
NSF-DMS 0409968. The first author was partially supported as a graduate stu- 
dent by the same grant. 

Department of Mathematics and Statistics, Dalhousie University, Nova Scotia, Canada 
B3H 3J5 

E-mail address: dmannaSmathstat.dal.ca 

Department of Mathematics, Tulane University, New Orleans, LA 70118 
E-mail address: vhmSmath . tulane . edu 



